Electron irradiation of Co, Ni, and P-doped 
BaFe2As2~type iron-based superconductors 



Cornelis Jacominus van der Beek , Sultan Demirdis, , Dorothee 
Colson 2 , Florence Rullier-Albenque 2 , Yanina Fasano 3 , Takasada 
Shibauchi 4 , Yuji Matsuda 4 , Shigeru Kasahara 45 , Piotr Gierlowski 6 , 
Marcin Konczykowski 1 

^aboratoire des Solides Irradies, CNRS UMR 7642 & CEA-DSM-IRAMIS, Ecole 
Polytechnique, F-91128 Palaiseau cedex, France 

2 Service de Physique de l'Etat Condense, CEA Saclay, CEA-DSM-IRAMIS, CNRS URA 2464, 
F-91191 Gif-sur-Yvette, France 

3 Laboratorio de Bajas Temperaturas, Centro Atomico Bariloche & Instituto Balseiro, Avenida 
Bustillo 9500, 8400 San Carlos de Bariloche, Argentina 

4 Department of Physics, Kyoto University, Sakyo-ku, Kyoto 606-8501, Japan 

5 Research Center for Low Temperature and Materials Sciences, Kyoto University,Sakyo-ku, 

Kyoto 606-8501, Japan 

institute of Physics of the Polish Academy of Sciences, 32-46 Aleja Lotnikow, 02-668 Warsaw, 
Poland 

E-mail: keesvanderbeekOpolytechnique . edu 

Abstract. High energy electron irradiation is used to controllably introduce atomic-scale point 
defects into single crystalline Ba(Fei_ ;E Co; C )2As2, Ba(Fei_ I Ni a; )2As2, and BaFe2(Asi_ a; P I )2. 
The appearance of the collective pinning contribution to the critical current density in 
BaFe2(Asi_ a; P I )2, and the magnitude of its enhancement in Ba(Fei_ a; Co a ,)2As2, conform 
with the hypothesis of quasi-particle scattering by Fe vacancies created by the irradiation. 
Whereas the insignificant modification of the temperature dependence of the superfluid density 
in Ba(Fei_ I Co a; )2As2 and Ba(Fei_ I Ni I )2As2 points to important native disorder present before 
the irradiation, the critical temperatures of these materials undergo a suppression equivalent to 
that observed in the much cleaner BaFe2(Asi_ a; P I )2. This lends credence to the hypothesis of 
line nodes of the order parameter (at finite fc z ) in the former two materials. 



1. Introduction 

The premise of s± superconductivity in the multiband iron-based superconductors [1, 2], 
with a sign-changing order parameter between the electron-like and hole-like Fermi-surface 
sheets [2, 3, 4], has raised the question of the effect of atomic-scale point-like disorder on 
superconductivity in these materials [5, 6]. In particular, interband scattering of quasiparticles 
by point-like disorder was suggested to be pair-breaking, and, therefore, to lead to a strong 
suppression of the critical temperature, the appearance of sub-gap states [7], and a peculiar T 2 
dependence of the London penetration depth at low temperature [8, 9, 10, 11, 12, 13, 14]. 

More recently, the question of the role of point-like disorder has become only more relevant 
[15]. This is because the anti-ferromagnetic coupling between Cooper pairs on different Fermi 
surface sheets [2, 3, 4], and the subtle changes of the Fermi surface with doping [3, 4], isovalent 



chemical substitution [17, 18, 19, 20, 21, 26, 27], and pressure [22, 23, 24, 25], allow for a wide 
variation of realizations of the order parameter [3, 4, 27]. In particular, the likelihood of line 
nodes of the order parameter in isovalently doped BaFe2(Asi_ a; P ;E ) 2 [32], LaFePO [29], and 
KFe2As2 [28, 30, 31] are indicative of a nodal s±- or a d-wave superconductivity [3, 4] in the 
first two and latter respective cases. 

Nodal gap structure not imposed by symmetry, but arising from e.g., the anisotropy of the 
Fermi surface and/or interband scattering [38], has also been evoked [16]. In the case of an s± 
state, the order parameter sign would not only change between bands but also within a given 
band [16]. One can also have an s-wave gap-structure in which nodes appear for non-zero k z 
[15, 39, 40], a situation that might be relevant for Ba(Fei_ x Coa;)2As2 [35], Ba(Fei_ :r Ni. r )2As2 
[12], Nai_5(Fei_ a; Co 3 ;)As [36], and overdoped Bai_. E K E Fe2As2 [37]. 

1.1. Disorder effect on exotic superconductivity in iron pnictides 

Given that the appearance of gap nodes in an s-wave superconductor does not imply 
symmetry breaking, indirect experiments yield, as yet, most insight into the order parameter 
structure. Such experiments include the low temperature behaviour of the penetration depth 
[8, 9, 10, 11, 12, 13, 14, 27, 29, 32] and the thermal conductivity [30, 32, 35, 37] as probes of the 
low-energy quasi-particle density of states (QPDOS), but also the sensitivity of the iron-based 
superconductors to quasi-particle scattering by homogeneous atomic-scale point-like disorder. 
This can be introduced by, e.g., chemical substitution or energetic particle irradiation. Namely, 
the occurrence of sign changes of the order parameter implies that even non-magnetic point-like 
scatterers may act as pair-breakers [5, 6, 7, 15]. For s± -superconductivity in particular, impurity 
scattering of quasi-particles between electron-like and hole-like bands (with opposite sign of the 
order parameter) would lead to a rapid suppression of the critical temperature T c [5, 6, 7], and 
open the possibility for non-sign changing s ++ superconductivity due to orbital fluctuations 
to manifest itself [5, 6, 41]. The vulnerability of the d-wave symmetry to both interband- 
and intraband impurity scattering would render this even more fragile than s±. On the other 
hand, a less pronounced role of interband compared to interband scattering would reduce the 
sensitivity of s± superconductivity to disorder [15, 41]. Finally, the presence of "accidental" 
nodes would result in a rather high initial sensitivity of superconductivity to disorder, with a 
slower suppression of the critical temperature as function of disorder strength once the interband 
scattering is sufficiently strong to wash out, or "lift" the nodes on all Fermi surface sheets [16]. 

1.2. Effect of native disorder 

Even if the problem at hand is characterized by a multitude of parameters, it is important to at 
least establish trends as controlled point-like disorder is added to iron-based superconductors. 
Early experiments have suggested a rather limited sensitivity of superconductivity in the iron- 
based materials to disorder. For example, the residual electronic term in the specific heat 
suggests a large low-temperature QPDOS [42, 43] away from optimal doping, that may be due 
to pair-breaking. The ubiquitous T 2 -behaviour of the low-temperature aft-plane penetration 
depth X a b(T) suggest strong pair-breaking [8, 9, 10, 11, 12, 13, 14]. Still, T c 's of the materials 
in question are high. 

In our earlier work [44, 45], we have drawn attention to the fact that vortex pinning in the 
mixed state of the iron-based superconductors is indicative of the importance of native disorder. 
While the critical current density maximum j c at low fields, and the subsequent j c ~ B~ x l 2 
drop-off are due to nm-scale heterogeneities of the superconducting properties, leading to a so- 
called "strong-pinning" contribution j s c to the critical current (see Fig. 1 and Refs. [47, 48]), the 
so-called "weak collective pinning" [49, 50] contribution j^ oU to the critical current density in the 
magnetic field range of several tenths of a T to several T is consistently interpreted in terms of 
quasi-particle scattering by the dopant atoms in charge-doped iron-based compounds. Estimates 



of scattering rates of the doping impurities, such as these are obtained from flux pinning, yield 
rather large values, which are at odds with a superconducting ground state that would be 
sensitive to point-like disorder should this be pair-breaking. Moreover, it may be noted that 
the occurrence of the quasi-particle scattering contribution to flux pinning consistently coincides 
with that of the T 2 behaviour of the penetration depth. Inversely, it is conspicuously absent in 
clean materials with a T-linear dependence of \ a bi such as BaFe2(Asi_ :E P x 

) 2 [32, 45, 58]. 

The relative insensitivity of superconductivity in the iron-based materials may be due to these 
materials not having a sign-changing order parameter, but may also be attributed to limited 
interband scattering by the dopant impurities. Given the premise of Coulomb scattering (with 
a small change in crystal momentum) by the charged dopants [45], this is not unreasonable. 
Estimates from pinning situate the native point defects rather in the Born limit (see table 1). 

1.3. Artificial disorder 

An obvious manner to test the above ideas is through the artificial introduction of different 
kinds of point-like defects, either by chemical doping [59, 60, 61], or by irradiation [62, 63]. 
An example is the substitution of magnetic ions, which suppresses T c rather more effectively 
than non-magnetic substitutions [59, 61]. A major drawback of chemical substitution in the 
iron-based compounds is, however, that this simultaneously leads to structural changes and/or 
doping. Energetic particle irradiation therefore seems preferable. Recent work by Tarantini 
et al. reported on T c -suppression by a-particle irradiation of NdFeAs(0,F), but was prone 
to criticism, notably in that it induces magnetic impurities and defect clustering. The 3 MeV 
proton irradiation experiments by Nakajima et al. on Ba(Fei_ ;E Co :E )2As2 reported monotonous 
T c -depression as the residual resistivity increases; the authors provided estimates for the critical 
pair-breaking parameter as T /2nT c ~ 4 — 7 and V /2ttT c ~ 1.5 — 2.5, as one goes from underdoped 
to overdoped [63]. However, apart from point defects, proton and neutron irradiation are known 
to induce, in situ, point defect cascades and clusters, which may play a different role than that 
of simple scatterers. 

Here, we compare the effect of 2.5 MeV electron irradiation on differently substituted 



Table 1. Contribution of dopant disorder to elastic scattering parameters of various iron 
pnictide superconductors, such as estimated from the weak collective pinning contribution to 
the (flux pinning) critical current density, j c ° u ■ Here, Hf is the Fermi wavevector, £o is the 
Bardeen-Cooper-Schrieffer coherence length, rid is the atomic point defect density, ot r is the 
transport scattering cross-section, 5o is the scattering phase angle, T = nd[7riV(0)] -1 sin 2 So is 
the scattering rate, where N(0) = rnkp/^TTh 2 is the normal-state DOS. 
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BaFe2As2- Such high-energy electron irradiation is known to produce homogeneously distributed 
vacancy-interstitial (Frenkel) pairs on all sublattices of a crystal, but, preferentially on that 
constituted by the lightest atomic species in the material - in the materials studied here, Fe. 
In the high-T c cuprates, electron-irradiation defects are known to be strong unitary scatterers, 
comparable to Zn substitution. There, they are responsible for the suppression of the critical 
temperature [64, 65, 66], and the appearance of a T 2 -temperature dependence of the penetration 
depth [66]. In semi-metals with extremely low carrier densities, such as Bi, the irradiation- 
induced defects are charged and lead to effective doping of the material [67]. The mobility of 
the created point defects at high temperature is a caveat, since it leads to partial annealing and 
clustering, albeit it to a lesser extent than that produced by proton or neutron irradiation. 
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Figure 1. 



Sin- 



gle crystalline 
Ba(Fe .925Co .o75)2As 2 : 
Hysteresis loop of the 
irreversible magnetiza- 
tion, at 5 K, before (o) 
and after (o) irradiation 
with 5.5 x 10 19 cm" 2 2.5 
MeV electrons. Closed 
arrows depict the di- 
rection in which the 
loop is traversed upon 
cycling the magnetic 
field. Dotted double 
arrows depict the width 
of the magnetization 
loop in the low-field, 
strong pinning regime, 
and in the higher field 
regime in which only 
the weak collective 
pinning contribution 
jcoll j g re i e vant. 



It turns out that electron irradiation leads to similar monotonic suppression of T c in 
Ba(Fei_a;Co x .)2As2, Ba(Fei_ a: Ni a ;)2As2, and BaFe2(Asi_ a; P :r )2- In Ba(Fei_ :r Co a ;)2As2, the 
suppression is weakest for the optimally doped material and strongest for underdoped crystals. 
In the range of intermediate electron fluences studied here, the resistance monotonically increases 
as function of fluence. In Ba(Fei_ x Coa;)2As 2 and Ba(Fei_ x Ni x )2As2, there is little to no change 
of the temperature dependence of the superfluid density. We have also measured the critical 
current density before and after the irradiation. In Ba(Fei_ a; Co 2: )2As2 one has a clear increase 
of the weak collective pinning contribution to j c , while in BaFe2(Asi_ x P x ) 2 the initially absent 
weak pinning contribution emerges after irradiation. This finding establishes that the field- 
independent contribution to the critical current density of iron-based superconductors is, 
indeed, due to atomic sized point pins. In both materials, this finding allows for an estimate of 
the density of produced defects per unit irradiation fluence. In Ba(Fei_ :r Ni x )2As2, j c is strongly 
suppressed by the irradiation. 




Figure 2. Transition from the superconducting to the normal state of a Ba(Feo.925Coo. 075)2^82 
single crystal before (a,c) and after irradiation with 2.1 Ccm~ 2 2.5 MeV electrons (b,d), as 
imaged by the differential magneto-optical (DMO) technique [47]. The images (a,b) show the 
progressive admission of an ac magnetic field of magnitude 1 Oe, applied perpendicularly to the 
sample surface, as the temperature is raised. In these flux density maps, areas of high luminous 
intensity / correspond to the value of the applied field, while dark areas correspond to zero 
field, (i.e. full screening). The dark rectangles in the upper left hand panels of (a) and (b) 
correspond to the sample outline, i.e., full Meissner expulsion of the magnetic field at the lowest 
temperature. The white squares indicate the positions at which the "transmittivity" (or "local 
ac susceptibility") defined as [I(T) - I(T < T C )]/[I(T > T c ) - I(T < T c )} depicted in the lower 
panels (c,d) was determined. 



2. Experimental details 

2.1. Single crystalline samples 

Single-crystals of Ba(Fei_ ;E Co a ;)2As2 [68] and Ba(Fei_ x Ni x )2As2 [69] were grown using the self- 
flux method. The high purity starting reagents Ba, FeAs, and CoAs were mixed in the molar 
ratio 1 : (4 — x) : x, loaded in alumina crucibles and then sealed in evacuated quartz tubes. 
For each doping level, chemical analysis by an electron probe was performed on several crystals 
yielding the dopant content within 0.5% absolute accuracy. For this work we studied six Co- 
doping levels and two Ni-doping levels. The optimally doped BaFe2(Asi_ x P x )2 single crystals, 
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Figure 3. Surface impedance Z s of Ba(Feo.925Coo.o75)2As2 crystals before and after various 
low-temperature irradiation runs with 2.5 MeV electrons, and subsequent annealing at 300 K. 
Data points and thin lines show the surface resistance R s and reactance X s respectively, for 
various electron doses, (b) Shift <5/(T) of the resonant frequency of the superconducting Nb 
cavity (operating at 26.0 GHz) as function of temperature, normalized with respect to the 
low-temperature extrapolated 5f(0), for Ba(Feo.925Coo.o75)2As2 crystals irradiated with various 
fluences of 2.5 MeV electrons. 

with x = 0.33 and x = 0.36, were also grown by the self flux method [70], and characterized 
using x-ray diffraction and energy dispersive x-ray spectroscopy. No impurity phases were found 
within the experimental limits of ~ %1. 

2.2. Electron irradiation 

The 2.5 MeV electron irradiation was performed at the SIRIUS Pelletron facility of the 
Laboratoire des Solides Irradies (LSI) at the Ecole Polytechnique in Palaiseau, France [71]. 
Beam currents varied between 10 and 20 /iA, depending on the irradiation run; the beam was 
swept over a 6 x 6 mm 2 area. Samples of different composition irradiated to the same dose were 
mounted together and irradiated simultaneously. In order to prevent in-situ defect migration, 
recombination, and clustering, the irradiations were carried out in a liquid H2 bath (T = 20 K). 
In-situ resistance measurements were carried out on the Ba(Fei_ x Co a ;)2As2 material in order to 
assess the increase of the residual resistivity as a function of defect density at low temperature, 
as well as the effect of room-temperature annealing. 

2. 3. Measurements of T c and the surface impedance 

T c -values before and after irradiation were determined using differential magneto-optical (DMO) 
imaging (see Fig. 2 and Ref. [47]). For the lower irradiation doses, in particular, the shift in T c 
after irradiation is smaller than the transition width. DMO then provides an objective means 
to quantify the effect of the electron irradiation on different portions of a crystal. The change in 
T c was further quantified using ex-situ measurements of the temperature dependent microwave 
surface impedance Z s = R s + iX s . These were performed at the LSI using a superconducting 
Nb cavity, operating at 26.0 GHz, and cooled to 5 K using a 0.5 W cryocooler cold head. From 



the temperature dependent resonance frequency / and unloaded quality factor Qq, the surface 
resistance R s and reactance X s were determined as 
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of the cavity, d the thickness of the platelet-like single crystal sample, and w its smaller width. 
The additive constant Xq was adjusted so that R s = X s in the normal state. 

2.4- Determination of the critical current density 

Measurements of the magnetization M of the Ba(Fei_ a; Co :r )2As2 and Ba(Fei_ ;E Ni :r )2As2 single 
crystals as function of the applied magnetic field H a were performed using a commercial 
SQUID magnetometer. From the width AM(H) of the magnetic hysteresis loops, the field- 
and temperature dependence of the sustainable current density j = |AM/u> was determined 
using the Bean model. The numeric prefactor | is estimated from calculations of E.H. Brandt for 
rectangular bars of similar aspect ratio as the measured crystals [72]. Data on BaFe2(Asi_ x P x )2 
were acquired by M. Konczykowski using the local Hall probe magnetometry technique [73, 74]. 
Briefly, the single crystalline superconducting crystal is centered on top of an array of 11 Hall 
sensors, of active area 3x3 /um 2 , and spaced by 20 (im, in such a manner that the short er crystal 
edge is parallel to the array, and the array spans the longer crystal edge. The Hall sensors are 
used to measure the spatial gradient dB/dx of the local induction component perpendicular to 
the crystal surface. The gradient dB/dx ( in G/^m ) is directly proportional to the sustainable 
screening current j, with dB/dx ~ gj for a crystal aspect ratio of ~ 0.1 [72]. Measurements of 
the flux density gradient were performed both as function of H a after zero field-cooling, and as 
function of temperature T, after field cooling, and the subsequent reduction of H a by twice the 
full penetration field. 

In what follows, we shall take the screening current at the lowest measurement temperature 
( 5 K ) to be representative of the critical current density. At higher temperature, the measured 
current density can be significantly lower than the pinning critical current density because of 
flux creep [58, 74]. 

3. Analysis of the critical current density 

The critical current density is further analyzed along the lines of Refs. [44, 45]. The hysteresis 
loops of all charge-doped iron-based superconductors, exemplified by Fig. 1, indicate the 
presence of two main pinning mechanisms in these materials. The central peak of the 
magnetization, at low applied fields, is due to the strong-pinning contribution ] s c to the critical 
current, while the constant contribution at higher fields comes from "weak collective pinning" 
by the dopant atoms. In a critical current density - versus magnetic field plot (see Fig. 6) 
the strong pinning contribution is responsible for the low-field plateau and the subsequent 
j c oc B^ 1 / 2 decrease, while the collective pinning contribution yields the high-field constant j c . 
It can be written in terms of superconducting parameters, defect density n^, and the elementary 
pinning force f p of a single defect, as 
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Figure 4. Ba(Fe .925Coo.o75)2As 2 : 
Resistance as function of dose of 2.5 
MeV electrons, measured at 21 K. The 
sample transits to the normal state 
after a dose of 1.2 Ccm -2 , after which 
the resistance increases at a rate of 
5R/R = O.WtCcm- 2 ]- 1 . The Inset 
shows the effect of annealing. The 
resistance of the crystal is represented 
as function of time. The protocol 
comprises initial measurements at 
21 K following irradiation with 4.6 
Ccm -2 2.5 MeV electrons. The 
crystal is then warmed to 300 K, 
cooled to 21 K, and heated once again 
to 300 K. After the first anneal at 
300 K, the resistance drops back to the 
value reached after low-temperature 
irradiation with only 1.4 Ccm -2 (white 
circle on the curve in the main panel). 



Here, jo is the depairing current density, £ a & and A Q & are the aib-plane coherence length and 
penetration depth respectively, e$ = ^q/Att^qX^ is the vortex line energy, E\ is the penetration 
depth anisotropy [51], and <l?o = h/2e is the flux quantum. The averaging (. . .) is performed 
over the vortex core. 

The magnitude of is compatible with scattering of the quasi-particles in the vortex 
cores as being at the origin of the weak collective pinning contribution in pristine iron-based 
superconductors, provided that the dopant atoms are the scattering defects. This gives rise to 
[52, 53] 

f p = OMPD)eo (^] (1°-) (4) 

which depends not only on the Gor'kov impurity parameter pp> = hvp/2irT c l = £o/h but 
also on the transport cross-section atr = (2ft /kp) sin 2 So (So is the scattering phase angle, and 
kp the Fermi wavevector). The analysis of the intermediate field critical current density of 
several common iron based superconductors allows one to estimate the scattering parameters 
compiled in Table 1. The weak collective pinning contribution is absent in isovalently doped 
BaFe 2 (Asi -nPi)2 for all doping levels x [see Ref. [58] and Fig. 9(a)]; surprisingly, it is present 
in Ba(Fe J _ :r Ru E )2As2 [54] and LiFeAs [73]. 



4. Results and Discussion 

4.1. Ba(Fei-. x Co x ) 2 As 2 

Fig. 4 shows the resistance of an optimally doped Ba(Feo.925Coo.o7s)2As2 single crystal, measured 
in- situ during exposure to the 10 fiA, 2.5 MeV beam, at T = 21 K. The sample transits 
to the normal state at a fluence of 1.2 Ccm -2 , after which the resistance increases at a rate 
of 5R/R = 0.14[Ccm -2 ] -1 . After a fluence of 4.6 Ccm -2 was reached, the electron beam 
was stopped, the sample heated to 300 K, and the resistance measured. The crystal was 
subsequently cooled down once again to 21 K, and the resistance measured again (Inset to 
Fig. 4). It appears that after the 300 K anneal, the resistance has dropped to the value first 
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Figure 5. Ba(Fei_ a: Co a .)2As2: evolution of T c (normalized by the initial critical temperature 
T c o) as function of (a) the fluence of 2.5 MeV electrons (b) the estimated normal state scattering 
rate V (normalized by T c q), for different doping levels x. The small data points show the T c -values 
after annealing at 300 K. The large green triangle shows the drop of T c following low-temperature 
irradiation of optimally doped Ba(Feo. 925C00. 075)2^82, the arrow shows the effect of annealing 
at 300 K. 



measured after irradiation with a fluence of 1.4 Ccm -2 . Therefore, annealing at 300 K leads to 
the annihilation and clustering of the point defects produced by the irradiation, with a resulting 
drop of the irradiation-induced resistance change of 65 %. Since all data presented below concern 
ex-situ experiments performed after heating the crystals to 300 K, the effect of annealing should 
be taken into account. 

The evolution of the reduced critical temperature T c /T c q with electron fluence (after 300 K 
annealing) is shown for various Co-substitution levels x in Fig. 5(a). Here, T c q is the critical 
temperature of the crystal before irradiation. It is to be stressed that, given the modest T c 
changes involved, and the typical dispersion of T c -values between similarly grown samples, even 
from the same batch, it is essential to compare the values of T c and T c $ on the same crystal 
before and after irradiation. For the same reason, a plot of the absolute T c versus dose would 
be meaningless in the regime of small doses. Fig. 5 also shows the effect of 300 K annealing, for 
one sample. It is seen that the change in T c drops by a factor three after annealing. 

Fig. 5(b) shows the same data, as function of the product of the induced change in the 
scattering rate 5T/2irT c and the electronic effective mass enhancement z = m*/m e . Here 
m e is the free electron mass, and m* the effective mass. z5T = e5p/ Rufn e (with e the 
electronic charge) is estimated from the change of the resistance due to the irradiation (after 
annealing), together with the published data for the resistivity p and the Hall coefficient .Re- 
published in Ref. [68]. In spite of the contributions of multiple bands to the conductivity of 
Ba(Fei_ x Co x )2As2, this approach can be justified given that the conductivity contribution from 
the electron sheets dominates the transport properties [68]. Typically, 5R/R ~ 0.05 [Ccm -2 ] -1 
for the optimally doped sample (cf Fig. 3). Note, however, that the resistivity change with 
electron dose depends on x, and that, in order to produce the "simulated" 5T values of Fig. 5b, 
we use data measured in-situ on crystals with different x. Estimates of the scattering rate from 
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Figure 6. Sustainable current density j of (a) optimally doped Ba(Fei_ x Coz)2As2 (x = 0.075) 
and (b) optimally doped Ba(Fei_ x Ni. r )2As2 ( x = 0.035) before (open symbols) and after 
irradiation at 21 K with 5.5 x 10 19 cm -2 2.5 MeV electrons (closed symbols). The measurements 
were performed at 5 K ( o ,•), 7.5 K (❖,<>), 10 K (□,■), 11 K (□,■), 13 K (A ), 16 K 

(V), and 17.5 K ( A , A ). 



critical current measurements (see below) yield qualitatively similar results as those shown in 
the Figure.. 

The surface impedance of optimally doped Ba(Fei_ a; Co a ;)2As2 before and after electron 
irradiation is depicted in Fig. 3. A monotonous increase of the surface resistance is observed after 
irradiation; even though in this series of experiments, the increase was not linear as function of 
fluence. In particular, the crystal irradiated with 2.1 Ccm -2 electrons did not conform to the 
general trend. The normalized frequency shift of the Nb cavity, which is proportional to the 
supernuid density n s oc A~ fe 2 of the superconducting crystals, showed a slight change to a less 
marked temperature dependence, indicative of disorder having been added to the crystal. The 
trend is similar to that observed by Hashimoto et al. in Bai_ x K a: Fe2As2 crystals with varying 
degrees of disorder [75], but is much less marked here. 

The low-temperature critical current density of single crystalline Ba(Feo.925Coo.o75)2As2 
increases by a constant, field- independent amount after electron irradiation (see Fig. 6a). This 
indicates that it is the weak collective pinning contribution that is enhanced, presumably 
by atomic-sized point defects introduced by the irradiation. If one assumes that the relevant 
defects are Fe vacancy-interstitial pairs, i.e. that atr ~ 7r Dp e ~ 2.6 A 2 , with Dp e = 0.92 A the 
ionic radius of Fe 2+ , the increase of the critical current density is consistent, through Eqs. (3) 
and (4), with 0.006 dpa Fe / Ccm~ 2 . The increase of the low-temperature critical current 
density and the concomitant decrease of T c after irradiation implies a steeper decrease of j c 
with temperature, and a crossing of j c -values of the pristine and the irradiated crystal at an 
intermediate temperature. 

4.2. Ba(Fe 1 ^ x Ni x ) 2 As 2 

Surface impedance measurements and DMO imaging on the irradiated Ba(Fei_ x Ni ;c )2As2 
crystals show a more rapid depression of T c with dose than in Ba(Fei_ a; Co :r )2As2. If anything, 




Figure 7. Temperature dependence of 
the flux density gradient dB/dx in a 
BaFe2(Aso.64Po.36)2 single crystal before and 
after irradiation with 0.8 x 10 19 electrons 
cnr 2 . Inset: Flux density profile across the 
crystal after zero-field cooling, application of 
the applied field ^H a = 100 mT, and subse- 
quent warming. 



Figure 8. Irreversibility field BHi rr (T) 
as function of temperature, for 
BaFe2(Aso.64Po.36)2 single crystals irra- 
diated with 0.8 x 10 19 (□), 1.6 x 10 19 (•), and 
4.5 x 10 19 electrons cm -2 (o, □, □). B irr (T) 
is determined as the locus of temperatures T 
above which the nonlinearity of the crystal 
current-voltage characteristic vanishes at 
given induction H a . T c is determined from 
the extrapolation of Hi rr (T) to zero field. 



T c versus dose is comparable to that found in underdoped Ba(Fei_ x Co x )2As2 (see Fig. 11). As in 
Ba(Fei_ x Co x )2As2, the temperature dependence of the superfluid density of Ba(Fei_ x Ni. c )2As2 
changes little or not at all, even at the largest irradiation dose of 5.5 x 10 19 electrons cm~ 2 . 

The critical current density of the Ba(Fei_ IE Ni x )2As2 crystals is strongly suppressed after 
the irradiation (see Fig. 6b). This is partially because the reduced measurement temperatures 
are higher in Ba(Fei_ x Ni x )2As2 than in its co-substituted counterpart, but mainly because the 
suppression of the prefactor jo oc n s e(A/hkF) in Eq. (3) outweighs the increase in rid (A is the 
superconducting gap amplitude). 

4.3. BaFe 2 (Asi- x P x ) 2 

The Inset to Fig. 7 shows an example of flux density gradients measured on a BaFe2(Aso.64Po.36)2 
single crystal irradiated with 0.9 x 10 19 electrons cm -2 , after zero field cooling and the application 
of an external field of 100 mT. The temperature dependence of the flux density gradient obtained 
after field-cooling in 400 mT, reduction of the applied field to the indicated target value, and 
subsequent warming is shown in the main panel of the Figure. The flux density gradient 
is directly proportional to the sustainable current density j. As in Ba(Fei_a;Co x .)2As2, the 
screening current in the irradiated crystal exceeds that of the pristine crystal at low temperature, 
but drops below it at higher temperature, due to the decrease of T c . 

The critical temperatures after electron-irradiation were determined as the extrapolation 
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Figure 9. (a) Hysteretic loops of the local flux density gradient versus local induction B, 
measured on the surface of a pristine BaFe2(Aso.67Po.33)2 single crystal, at various temperatures 
(indicated). The astroid-shaped hysteresis loops are determined by the sole strong-pinning 
contribution to the critical current, (b) ibid, measured on the surface of a BaFe2(Aso.7Po.3)2 
single crystal irradiated at 23 K with 1.8 x 10 19 electrons cm -2 , at the indicated temperatures. 
The opening of the loops at higher flux densities reveal the emergence of a weak collective 
pinning contribution by the atomic-sized point defects introduced by the irradiation. 

to zero of the "irreversibility field" Hi rr (T) above which the nonlinearity of the current-voltage 
characteristic - and therefore the critical current density - vanishes. Details on the determination 
of Hi rr , can be found in, e.g., Ref. [74]. In all cases, the Hi rr (T) values depend very little on the 
frequency of the ac field used to investigate the screening by the superconducting sample; the 
corresponding Ti rr {H) lies very close to the temperatures at which the dc screening current 
vanishes in Fig. 7. The T c -values resulting from the extrapolation of Hi rr {T) to zero are 
gathered in Fig. 11. The T c -depression as function of electron fluence is comparable to that in 
Ba(Fei_ a; Ni 2: )2As2. Note that the resistivity increase of the irradiated BaFe2(Asi_ x P x )2 crystals 
corresponds to 6R/R ~ 0.16 [Ccm" 2 ]" 1 , three times higher than for Ba(Fei_ x Co x )2As2. 

Figure 9 shows the hysteretic loops of the local flux density gradient on the surface of a pristine 
(a) and electron-irradiated BaFe2(Asi_ x P 2: )2 single crystal (b), for different temperatures. The 
pristine single crystal shows astroid-shaped magnetic hysteresis loops that are indicative of 
strong flux pinning by extended point-like defects only. After electron irradiation, the hysteresis 
loops open up due to the appearance of the field-independent collective pinning contribution 
from atomic-sized point defects. We interpret the appearance of a non-zero j™ as being 
due to vortex pinning by vacancy-interstitial pairs introduced by the irradiation. Moreover, 
jcoii m onotonically increases (at low temperature) as function of irradiation fluence (Fig. 10), 
which allows for a direct comparison with the theory for quasi-particle scattering mediated 
collective vortex pinning. Under the assumption that the only unknown parameter, at r , is 
again approximated by the ionic cross-section of a Fe-vacancy, a fit to Eq. (3) yields a defect 
density of 0.0035 dpa / Ccm -2 . This number is taken as more precise than that obtained for 
Ba(Fei_ x Co x )2As2 in subsection 4.1, in which only two points were available. Note that the 
hypothesis of simple voids (non-scattering point defects) would necessitate an unphysical 1 dpa 
/ Ccm~ 2 to explain the magnitude of the critical current density change. 



4-4- Discussion; Comparison between differently substituted materials 

The modification of the critical current density of Ba(Fei_ x Co x )2As2 and BaFe2(Asi_a;P x )2 by 
electron irradiation demonstrates that the field-independent contribution to j c apparent in sub-T 
to T-fields is indeed due to collective pinning by atomic-sized point defects. The very reasonable 
number obtained for defect generation by the electron beam, 0.35% dpa / Ccm~ 2 , indicates 
that it is, in a first approximation, safe to assume that the most relevant produced defects 
are Fe vacancies, and that these vacancies are responsible for quasi-particle scattering. This 
interpretation is strengthened by the fact that scattering rates in Ba(Fei_ a: Co ;E )2As2 estimated 
from the change of the resistivity are very much comparable to T = n^/lirN^O)]/ sin 2 5q estimated 
from the defect density. A comparison with Table 1, shows that the scattering cross-section 
of Fe vacancies should be comparable to that of Co and Ru impurities. Indeed, electron- 
irradiation of Ba(Fei_ x Coa;)2As2 and Ba(Fei_ :r Ni ;E )2As2 does not lead to appreciable changes 
of the temperature dependence of the superfluid density of those materials, which suggests that 
the effect of the dopant impurities overwhelms that of the defects added by the irradiation. 

Irrespective of the precise changes in resistivity incurred by the irradiation, a striking 
feature of our results is the similar suppression of T c by the electron irradiation of underdoped 
Ba(Fei_a;Co x )2As2, Ba(Fei_ :r Ni. r )2As2 and BaFe2(Asi_ x .Pa;)2 for a given dose. The T c - 
suppression in optimally doped Ba(Fei_ x Co 1 .)2As2 is a factor 2 less. Given the clear evidence 
for line nodes in BaFe2(Asi_ x P a .)2, this poses the question possible nodal structures of the order 
parameter in Ba(Fei_ a; Ni a ;)2As2, and, possibly, Ba(Fei_ x Co x )2As2 as well. In the latter two 
materials, it was proposed that a nodal line might exist on the a-(hole-like) sheet, at finite k z 
[35, 12]. 

It is tempting to compare our results to that obtained by other types of irradiation, as well 
as by doping. It turns out that the suppression of T c for optimally doped Ba(Fei_a;Co x )2As2 
is remarkably comparable to that reported for Co-doping of Bai_ a; K E Fe2As2 in Ref. [76], 
and much less than that obtained by Co-doping of KFe2As2 studied by the same authors. 



Figure 10. 

BaFe 2 (Asi_ :r P :E )2 
(with x = 0.33, 
x = 0.36): Electron- 
fluence dependence 
of the weak collective 
pinning contribution 
j%° u to the critical 
current density. The 
drawn line denotes a 
fit to Eqs. (3,4). As- 
suming that the point 
defects most relevant 
for quasi-particle scat- 
tering are Fe vacancies 
induced by the irradi- 
ation, the parameter 
value 5.1 x 10 s would 
correspond to 0.0035 
dpa / Ccm -2 . 
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Figure 11. (a) Dose-dependence of the T c -change of Ba(Fei_ a; Co a ;)2As2, Ba(Fei_ x Ni a; )2As2, 
and BaFe2(Asi_ s P x )2, after irradiation with 2.5 MeV electrons and annealing at 300 K. 
The upper scale shows the density of point defects (presumably Fe vacancies) added by the 
irradiation, such as determined from the dose-dependence of the critical current density of 
BaFe2(Asi_a;P x )2 (see Fig. 10). (b) Dose dependence of T c /T c q, where T c q is the critical 
temperature of the pristine crystal. 



The decrease of T c observed after electron irradiation is also comparable, as far as optimally 
doped Ba(Fei_ x Co x )2As2 is concerned, to the Co-doping experiments on Bai_ 1 .K r Fe2As2 by 
Li et al, while the T c suppression in underdoped Ba(Fei_ ;c Co a ;)2As2, Ba(Fei_ rr NI r )2As2, and 
BaFe2(Asi_a;P x )2 is similar to the Cu and Zn doping by those authors. Even if the comparison 
with chemical doping is tenuous, the observed trend clearly indicates a much weaker sensitivity 
of the materials studied here than what is expected for the scenario of s± superconductivity with 
strong interband scattering, and sets the three materials even further from the <i-wave scenario. 

A comparison with the previous irradiation studies reveals that T c -suppression after electron 
irradiation is somewhat weaker than the results obtained by Nakajima by 3 MeV proton 
irradiation of optimally doped Ba(Fei_ a; Co :r )2As2 [63], and much weaker than in the a-particle 
irradiation of NdFeAs(0,F) by Tarantini et al. [62]. 

5. Summary and Conclusions 

Ba(Fei_ x Co x )2As2 crystals of different dopng levels x, as well as Ba(Fei_ x Ni x )2As2 and 
BaFe2(Asi_ :r P 2: )2 crystals, have been subjected to high energy (2.5 MeV) electron irradiation 
at low temperature (21 K). In spite of substantial annealing and clustering of the defects, 
the enhancement of the collective pinning contribution to the critical current density due 
to quasi-particle scattering in the vortex cores in Ba(Fei_ x Co x )2As2, and its appearance in 
BaFe 2 (Asi_ :r P x )2 in which it is absent in the pristine material, demonstrate that the irradiation 
produces atomic-sized point defects. The increase in magnitude of the low-temperature critical 
current density is consistent with the pinning by Fe vacancies, created by the irradiation. In 
Ba(Fei_a;Co x )2As2 and Ba(Fei_ x Ni. c )2As2, the introduction of these supplementary point defects 
does not lead to a significant modification of the temperature dependence of the superfluid 
density, which suggests that the defects are overwhelmed by pre-existing disorder (notably, the 
presence of the dopant atoms). Surprisingly, the critical temperature is similarly suppressed in 
all three materials. This, in spite of the fact that the order parameter in BaFe2(Asi_ :r P 2: )2 is 



thought to have line nodes, while this possibility is much less certain in the other two materials. 
The results therefore lend credence to evidence for line nodes obtained from c-axis penetration 
depth [12] and thermal conductivity measurements [35]. 
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